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GEOGRAPHICAL  COORDINATE  COMPUTATIONS 
Part  I 

General  Considerations 
By  W.  R.  Tobler 
Technical  Report  Number  2 
December,  1964 

abstract 


Part  I  provides  a  discussion  of  the  usefulness  of  coordinate  models 
for  studies  of  geographically  distributed  phenomena  with  comments  on 
specific  coordinate  systems  and  their  relevance  for  the  analysis  and  in¬ 
ventorying  of  geographical  information.  Appendices  include  equations  for 
conversion  from  the  Public  Land  Survey  system  into  latitude  and  longitude 
and  to  rectangular  map  projection  coordinates.  Part  II  considers  map 
projections  in  greater  detail,  including  estimates  of  the  errors  introduced 
by  the  substitution  of  map  projection  coordinates  for  spherical  coordinates 
Statistical  computations  of  finite  distortion  are  related  to  Tissot  s  In- 
dicatrix  as  a  general  contribution  to  the  analysis  of  map  projections. 
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INTRODUCTION 


In  recent  years  there  has  been  a  rapid  increase  in  the  use  of  formal 
mathematical  and  statistical  methods  for  the  analysis  of  terrestrial  dis¬ 
tributions.  Such  procedures  have  been  found  to  be  of  considerable  assist¬ 
ance  in  fields  such  as  city  and  regional  planning,  demography,  ecology, 
geography,  geology,  and  regional  science.  The  present  study  is  concerned 
with  only  one  of  the  several  mathematical  strategies  which  have  been 
utilized  for  such  analyses;  the  "coordinate  model".  This  term  is  taken 
to  include  that  class  of  studies  which  specifically  refers  to  the  location 
of  observational  phenomena  by  some  system  of  coordinates. 

As  an  example,  a  technique  associated  with  contemporary  theories  in 
geology  consists  of  estimating  the  departures  of  empirical  geological 
observations  from  a  "regional  trend".  Here  one  has  a  collection  of 
numerical  observations  (z^)  at  specific  terrestrial  locations  (x^,  y^), 
i  -  1,2,  ...  ,  n.  The  procedure  begins  by  estimating  a  specific  portion 
of  the  locational  trend  of  the  observations  as  a  least  squares  equation 
z  —  f(x,  y) ,  The  trend  is  then  subtracted  from  the  observations  to  obtain 
the  residual,  which  is  subsequently  interpreted  in  terms  of  geological 
theory.  In  this  instance  the  recording  of  locations  in  some  system  of 
coordinates  is  an  essential  prerequisite  for  the  analysis.  There  are 
many  other  such  examples.  The  coordinate  model  is  widely  used,  and  is 
enjoying  increasing  popularity  since  the  advent  of  electronic  computers, 
which  permit  facile  manipulation  of  the  metricized  locations.  The  re¬ 
searcher  now  also  has  available  to  him  a  rapidly  increasing  amount  of 
information  recorded  in  terms  of  coordinates;  the  U.S.  Bureau  of  the 
Census,  for  example,  now  provides  population  statistics  in  terms  of 
latitude  and  longitude  coordinates.  It  is  difficult  to  over-estimate  the 
usefulness  of  this  manner  of  recording  information  since  a  large  number 
of  the  analytical  methods  designed  for  the  analysis  of  distributions 
assume  the  existence  of  a  system  of  coordinates. 

As  a  system  of  locational  labels  the  specific  coordinates  employed 
for  the  recording  of  observations  are  not  of  direct  or  inherent  interest, 
but  rather  only  what  they  enable  one  to  deduce  regarding  interrelations 
among  the  phenomena  observed.  In  this  sense  the  particular  coordinate 
system  utilized  is  irrelevant.  On  the  other  hand,  computations  may  often 
be  simplified  by  the  choice  of  a  convenient  coordinate  notation.  From  a 
scientific  point  of  view  descriptions  of  phenomena  and  their  interrelations 
are  often  simplified  by  appropriate  formalizations  involving  a  "natural" 
coordinate  system  for  that  phenomena;  the  use  of  geomagnetic  coordinates 
in  the  study  of  terrestrial  magnetism,  for  example.  The  present  study, 
hOT7ever,  considers  only  systems  v;hich  appear  to  be  of  practical  utility 
for  the  large  scale  recording  of  terrestrial  observations,  with  some 
emphasis  on  systems  available  in  the  United  States. 

The  actual  surface  of  the  earth  can  be  referred  to  as  the  topographic 
surface.  This  bumpy  two-dimensional  surface  is  difficult  to  describe  in 
all  its  detail.  Theoretically  it  is  possible  to  introduce  a  system  of 
coordinates  on  this  surface  such  that  ground  distances,  etc.,  between  all 
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points  can  be  calculated.  In  practice  this  is  not  attempted.  As  an 
approximation  to  the  topographic  surface  the  geodesist  utilizes  a  surface 
of  constant  gravitational  potential,  the  geoid.  This  bumpy  but  rather 
smooth  surface  is  still  too  complicated  for  practical  computations.  A 
further  simplification  is  made  by  assuming  the  earth  to  have  the  shape 
of  an  ellipsoid,  generally  an  ellipsoid  of  revolution.  Geodetic  coordinates 
are  then  defined  for  positions  on  this  ellipsoid.  An  even  simpler  model 
of  the  topographic  surface  is  to  consider  the  earth  to  be  a  perfect  sphere. 
One  can  continue  thusly,  finally  arriving  at  the  assumption  that  the 
earth  is  a  flat  plane.  Each  of  these  assumed  models  of  the  earth  has  its 
advantages  and  disadvantages  since  realism  may  lead  to  extreme  cumbrousness. 
In  practical  terms,  the  following  (somewhat  contradictory)  criteria 
seem  appropriate: 

a)  The  coordinates  should  permit  accurate  and  economical  formulae 
for  computation.  The  highest  level  of  precision  available  today  can  be 
achieved  only  through  the  use  of  geodetic  formulae.  These  formulae  are 
fairly  complicated.  Computational  simplicity  can  be  obtained,  with  a 
consequent  reduction  in  accuracy,  by  employing  spherical  formulae. 

Further  computational  simplicity  can  be  achieved  by  the  use  of  plane 
coordinates  based  on  an  appropriate  map  projection,  again  with  some 
loss  of  accuracy. 

Computational  simplicity  is  important  for  two  reasons.  The  cost 
and  time  required  for  computation,  particularly  when  large  amounts  of 
information  are  to  be  processed,  can  be  reduced  by  significant  amounts 
through  the  use  of  simplified  formulae.  In  addition,  the  number  of 
persons  and  agencies  who  can  effectively  make  use  of  the  information  is 
significantly  increased  by  the  use  of  simplified  formulae. 

It  is  more  difficult  to  discuss  accuracy,  upon  which  the  level  of 
computational  simplificty  depends,  since  the  degree  of  precision  required 
in  particular  studies  varies  considerably.  There  is  no  reason  to  employ 
a  method  which  results  in  accuracies  greater  than  required  or  greater 
than  those  with  which  the  information  was  recorded.  An  objective  of 
this  study  has  been  to  estimate  the  accuracy  obtainable  by  employing 
several  alternate  methods.  This  allows  the  individual  researcher  to 
choose  the  simplest  computational  method  which  yields  the  requisite 
level  of  accuracy. 

b)  A  rapid  and  accurate  method  of  determining  the  coordinates  of  a 
position  should  be  available.  In  general  a  system  of  coordinates  which 
requires  a  carefully  executed  geodetic  survey  can  be  considered  highly 
accurate,  but  relatively  slow.  A  system  which  enables  one  to  read 
coordinates  from  an  aerial  photograph  or  map  (either  manually,  mechanically, 
or  electronically)  is  more  rapid  but  the  accuracy  is  dependent  on  the 

map  scale.  The  convenience  of  this  method  is  also  dependent  on  the 
availability  of  maps  or  photographs  to  which  the  coordinate  system  has 
been  affixed.  Between  the  extremes  of  geodetic  surveying  and  map  scaling 
are  a  number  of  intermediate  systems,  including  automatic  navigation 
devices  which  permit  virtually  instantaneous,  Jui  situ  position  determination, 
with  fair  accuracy.  Emphasis  in  this  study  is  on  map  scaling  procedures. 
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c)  The  coordinates  should  be  widely  available  and  should  be  equally 
convenient  for  use  at  a  local,  national,  and  international  level.  This 
objective  arises  since  most  types  of  information  collected  at  a  national 
level  are  used  both  nationally  and  locally.  Census  records  provide  a 
good  example.  National  use  of  local  records  also  is  increasing.  The 
accuracy  requirements  at  these  two  levels  generally  differ,  however. 

At  the  local  level  an  accuracy  of  fifty  feet  may  be  insufficient,  whereas 
at  the  national  level  an  accuracy  of  five  miles  may  suffice. 

terrf-strtat.  coordinate  systems 

There  are  many  locational  coordinate  systems  in  use  throughout  the 
world.  Emphasis  in  the  current  discussion  is  on  systems  available  in  the 
United  States. 

Geodetic  Coordinates; 

Geodetic  latitude  and  longitude  provide  the  traditional  method  of 
identifying  locations  on  the  surface  of  the  earth.  The  earth  is  assumed 
to  be  an  oblate  spheroid  and  the  geodetic  coordinates  are  based  on  actual 
measurement  (triangulation)  between  sets  of  locations  on  the  topographic^ 
surface.  These  values  are  then  adjusted  to  fit  an  ellipsoid  representative 
of  the  region  in  question.  Different  ellipsoids  are  employed  for  the 
several  continents  of  the  world,  with  an  International  Ellipsoid  in  use 
for  world-wide  computations.  Geodetic  coordinates,  based  on  the  Clarke 
Ellipsoid  of  1866  (1927  adjustment),  are  indicated  on  maps  published  by 
the  U.S.  Geological  Survey  and  by  the  U.S.  Coast  and  Geodetic  Survey. 

Latitude  and  longitude  scaled  from  such  maps  are  geodetic  coordinates, 
but  such  scaling  will  not  yield  the  same  accuracy  as  when  the  positions 
are  established  in  the  field  by  an  expensive  first  order  geodetic  survey. 

Computations  employing  geodetic  coordinates  usually  take  into  account 
the  ellipsoidal  shape  assumed  for  the  earth.  The  relevant  formulae  are 
fairly  complicated.  For  precise  geodetic  work  it  is  necessary  to  carry 
approximately  fifteen  significant  digits.  Experience  with  a  digital 
computer,  however,  indicates  that,  once  programmed,  the  ellipsoidal 
formulae  do  not  require  appreciably  more  effort  than  the  simpler  spherical 
formulae.  A  floating  point  program  with  seven  significant  digits  (as  employed) 
for  this  study)  yielded  values  which  differed  less  than  100  meters  from  more 
precise  values  over  a  range  of  6000  kilometers. 

Assumptions  required  to  apply  geodetic  computations  to  the  surface  of 
the  earth  are  that  (a)  the  geodetic  latitudes  and  longitudes  are  known 
without  error,  (b)  the  ellipsoid  chosen  is  representative  of  the  region 
in  question,  and  (c)  the  points  involved  lie  on  the  surface  of  the  ellipsoid 
(roughly,  at  sea  level).  On  the  other  hand,  this  is  the  most  accurate 
system  available.  The  actual  proportional  error  in  distance,  based  on 
misclosures  of  the  U.S.  continental  triangulation  network,  appears  to  be 
on  the  order  of 

_ 1 _  d1/3 

20000 

where  D  is  the  computed  distance  in  miles  on  the  Clarke  Ellipsoid  of 
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1866  (1927  adjustment).  The  differences  between  the  several  ellipsoids 
in  use  throughout  the  world  are  small;  on  the  order  of  three  kilometers 
per  6000  miles.  Connections  of  this  length  on  one  ellipsoidal  datum 
are  rare  and  the  figure  given  does  not  take  into  account  the  fact  that 
the  relation  between  the  several  datums  in  actual  use  are  not  yet  known 
in  detail;  in  other  words,  distances  between  positions  whose  geodetic 
coordinates  are  referred  to  different  ellipsoids  may  be  in  error  by  a 
larger  figure. 

Astronomic  Coordinates; 


Astronomic  latitude  and  longitude  are  based  on  celestial  observations 
and  may  depart  from  geodetic  coordinates  by  as  much  as  two  kilometers 
at  any  point,  due  to  departure  of  the  geoid  from  the  ellipsoid.  Astronomical 
observations  are  usually  available  only  for  isolated  points,  and  will  not 
be  considered  in  this  report. 

The  U.S.  Public  Land  Survey 

The  Public  Land  Survey  system  is  based  on  a  set  of  six  mile  squares 
numbered  as  townships  north  and  south  of  a  base  parallel,  and  as  ranges 
east  and  west  of  a  base  meridian.  These  six  mile  squares  are  then  sub- 
divided  into  36  sections,  each  one  mile  square,  and  numbered  in  serpentine 
fashion.  Each  section  can  be  further  subdivided  into  quarter-sections, 
each  one  sixteenth  of  a  mile  in  area.  Several  systems  similar  to  the 
Public  Land  Survey  exist  in  various  parts  of  the  United  States;  these 
are  not  considered  here. 

Strictly  speaking,  the  Public  Land  Survey  is  an  areal  identification 
scheme  and  not  a  metrical  coordinate  system,  though  it  is  often  regarded 
as  such.  As  a  partitioning  of  areas  the  system  does  not  differ  from 
county  or  census  units,  except  that  the  elemental  areas  are  roughly  of 
equal  size  and  are  labeled  in  a  more  convenient  fashion.  The  system  is 
not  complete,  in  the  sense  that  it  is  defined  only  for  certain  portions 
of  the  western  United  States.  In  these  areas  large  amounts  of  information 
have  been  (and  continue  to  be)  collected  and  recorded  in  terms  of  Sections, 
Townships,  and  Ranges.  These  collections  of  information  provide  a  valuable 
source  of  raw  data  for  research  workers.  The  Public  Land  Survey,  however, 
was  not  designed  for  the  analytical  manipulations  usually  required  in 
research  work.  For  example,  statistical  analyses  of  spatial  distributions 
may  require  calculation  of  the  average  location  (and  its  variance,  and 
so  on)  of  phenomena.  For  such  computations  the  distances  between  observed 
locations  may  be  needed.  The  distance  between  the  SW^,  Sec.  25,  T5S,  R7# 
Willamette  Meridian,  and  the  NE^,  Sec.  2,  T6N,  R8E,  Black  Hills  Merikan,’ 
is  not  immediately  apparent,  nor  is  there  any  simple  formula  which  can  be 
employed  to  obtain  this  difference.  Observations  recorded  in  the  Public 
Land  Survey  System,  however,  can  be  convered  to  a  coordinate  system  having 
the  requisite  metrical  properties.  This  is  because  the  Public  Land  Survey 
has  many  of  the  topological  ordering  properties  of  a  coordinate  system. 

The  most  direct  and  convenient  conversion  is  to  latitude  and  longitude. 

This  can  be  effected  in  several  ways.  The  system  of  Townships  and  Ranges 
is  shown  on  U.S.  Geological  Survey  topographic  maps  and  approximate 
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coordinates  could  be  scaled  from  there  maps,  A  more  convenient 
procedure  is  to  attempt  a  direct  calculation.  The  equations  for  such 
a  conversion  are  given  in  the  Appendix,  along  with  an  estimate  of  their 
validity.  The  errors  are  fairly  small  so  that  they  might  be  of  little 
consequence  when  working  with  observations  from  the  entire  United 
States,  The  urban  researcher  working  within  one  city,  on  the  other 
hand,  might  find  these  errors  intolerable. 

The  GEOREF  and  Marsden  Squares  Systems 

The  GEOREF  System  is  used  by  the  U,S,  Air  Force  to  identify 
locations.  It  is  a  modification  of  latitude  and  longitude  in  which 
letters  are  substituted  for  the  numerical  values.  Every  combination  of 
letters  is  taken  to  represent  a  quadrilateral  bounded  by  latitude 
and  longitude.  In  this  sense  the  system  is  a  partitioning  of  area 
rather  than  a  true  coordinate  system.  The  same  results  can  be  achieved 
by  using  latitude  and  longitude  with  a  convention  regarding  the  quadrant 
in  which  the  quadrilateral  of  area  lies.  The  system  has  certain  advantages 
in  applications  which  require  error- free  rapid  verbal  communications  (e,g, 
radio).  The  system  is  shown  on  maps  published  by  the  U,S,  Air  Force, 

The  Marsden  Squares  system  employed  by  the  National  Oceanographic 
Data  Center  is  similar  to  the  GEOREF  System  in  that  a  numbering  of 
]_atitude  and  longitude  quadrilaterals  is  substituted  for  the  geodetic 
coordinates.  There  are  many  other  such  systems  available,  including 
the  World  Aeronautical  Chart  designations  and  the  International  Millionth 
Map  of  the  World  system.  The  advantage  of  these  systems  is  largely  one 
of  bookkeeping.  Such  systems  are  not  further  considered  in  this  report. 

THE  SPHERICAL  ASSUMPTION 

The  various  computations  are  simplified  if  it  is  assigned  that  the 
earth  is  a  sphere.  The  results  of  such  computations  do  not  differ  by 
large  amounts  from  the  corresponding  ellipsoidal  values  -  the  polar 
flattening  of  the  earth,  after  all,  is  quite  small.  On  the  basis  of  a 
number  of  computations  it  appears  that  a  reasonable  and  convenient  rule 
of  thumb  is  that  the  flattening  of  the  earth  can  be  taken  as  an  approxjjna_^ 
upper  bound  on  the  percentage  error  of  measures  calculated  on  a  spherical 
as  compared  to  an  ellipsoidal  assumption.  This  is  about  one  part  in  300. 

An  even  safer  estimate  is  that  the  error  will  be  less  than  one  percent. 

For  some  purposes  this  is  intolerably  large,  but  for  the  majority  of 
requirements  it  is  far  more  accurate  than  are  the  data  or  theories  now 
available.  Detailed  numerical  differences  between  an  ellipsoid  and 
sphere  for  distances  and  angles  also  have  recently  been  published. 
Computation  of  the  differences  for  a  random  sample  of  200  pairs  of  points 
within  the  continental  United  States  resulted  in  the  following  values; 


Angular  differences  (degrees 


Distance  Differences  (miles) 
Average:  0.046 
Standard  Deviation:  1.896 
Minimum:  -3.788 

Maximum:  4.871 


Average:  0.006 
Standard  Deviation:  0.083 


Minimum:  -0.150 


Maximum:  0.159 


A  second  sample  might  yield  somewhat  different  results,  but  the  sample  is 
probably  representative  for  the  country  as  a  whole.  As  expected,  the 
differences  depend  on  both  the  distance  and  on  the  direction  of  the  point 
pairs.  A  comparison  of  surface  areas  is  given  in  the  accompanying  table. 

In  performing  these  computations  it  has  been  assumed  that  the  earth 
is  a  sphere  whose  radius  is  equal  to  the  equatorial  radius  of  the  Clarke 
Ellipsoid  of  1866,  and  that  the  geodetic  latitude  and  longitude  can, 
without  modification  of  their  numerical  values,  be  considered  to  be 
spherical  coordinates.  These  assumptions  have  the  advantage  of  extreme 
simplicity.  A  slight  improvement  in  accuracy  can  be  obtained  if  they 
are  not  retained.  For  example,  the  spherical  radius  employed  might  be 
the  average  radius  of  terrestrial  ellipsoid  in  the  vicinity  of  the  area 
of  interest,  rather  than  the  equatorial  radius.  It  can  be  proven  that 
the  average  radius  at  any  latitude  is  the  geometric  mean  of  the  radii  of 
curvature  along  the  meridian  and  normal  to  the  meridian.  This  average 
radius  is  given  in  the  accompanying  tables.  Conversion  of  geodetic 
latitude  to  spherical  latitude  can  also  be  accomplished  in  a  large  number 
of  ways.  Four  of  the  simpler  methods  are  illustrated  in  the  figure. 
Mathematical  treatments  can  be  found  in  works  on  geodesy  and  map  projections. 

THE  PLANE  ASSUMPTION 

It  often  is  convenient  to  employ  plane  coordinates  for  the  inventorying 
of  analysis  of  terrestrially  distributed  phenomena.  In  particular,  many 
of  the  numerous  statistical  and  analytical  methods  which  have  been  de¬ 
vised  for  the  analysis  of  two  dimensional  distributions  assume  the  exis¬ 
tence  of  a  system  of  Cartesian  coordinates.  As  a  very  simple  example, 
suppose  that  an  objective  is  to  compute  the  average  location  and  the 
locational  variance  of  a  set  of  discrete  phenomena  on  the  surface  of  a 
sphere.  One  can  proceed  in  several  ways; 

a)  Record  the  observations  in  latitude  and  longitude  and  then  perform 
the  calculations  using  the  spherical  formulae  for  average  and 
variance. 

b)  Plot  the  distribution  on  a  map,  assign  arbitrary  rectangular 
coordinates  to  the  map,  record  the  observations  in  these 
coordinates,  and  then  perform  the  calculations  using  the  plane 
formulae  for  the  average  and  variance. 


COMPARISON  OF  AREAS  FOR  A  ONE  DEGREE  ZONE 
OF  LONGITUDE  WITHIN  THE  UNITED  STATES 


(Values  in  square  miles,  rounded  to  the  nearest  square  mile) 


Latitude 


Ellipsoidal  Area  Spherical  Area* 


26N 

to 

27N 

4265 

4282 

27N 

to 

28N 

4228 

4244 

28N 

to 

29N 

4189 

4205 

29N 

to 

30N 

4150 

4164 

30N 

to 

31N 

4109 

4123 

31N 

to 

32N 

4067 

4080 

32N 

to 

33N 

4024 

4035 

33N 

to 

34N 

3979 

3990 

34N 

to 

35N 

3934 

3943 

35N 

to 

36N 

3887 

3895 

36N 

to 

37  N 

3839 

3846 

37N 

to 

38N 

37  89 

3796 

38N 

to 

39  N 

37  39 

3744 

39N 

to 

40N 

3687 

3692 

40N 

to 

41N 

3634 

3638 

4  IN 

to 

42N 

3581 

3583 

42N 

to 

43N 

3526 

3528 

43N 

to 

44N 

3469 

3471 

44N 

to 

45N 

3412 

3413 

45N 

to 

46N 

3354 

3354 

46N 

to 

47N 

3295 

3294 

47N 

to 

48N 

3234 

3232 

*Radius  equal  to  equatorial  radius  of  Clarke  ellipsoid  of  1866 


^0'0«o^tt<#^tf^^-<cw^•^cDCD^>o^'Or><#(Maomu^^-oo^^llr^f<>u^oo^^-'0^^ooo 

Or*-0'0’«tO^-'0<^4^^^4'OCO<#<M(<>CO<MOO'O^MJBOrgO^••^<MO'm^-^0^•M^O 

(M^OOO'0<MflO^«C^<0(*^II^U^U^«D^-^>00'^^*fM^OOO'OCDOO'fMf*^^-^0 

^-^»u^^4ao<^<o«DO^<^eo<0(^JCDr<^<AO*>o^^OO''0(<>^tf^O'•fcDrMu^ao^^^u^^- 
<0  'O'O  'O'O'O  'C  'Ou^tAu^u^tf^u^u^u^<^'^<#<#*t'r<^r<^^^ 


r*^i0O'0'^or^'0Oor^o®0'O®O(M0^h*<MO'0«0h-r^rM^’^^^Of^0'^0'fM 

«D<0^^®0'0^><^J<^^>^>U^(M^><0^®<OU^(*>®®f*>^^<^ifMO'0®'0®U^^-^0 

U^mU^U^tAU^U^'0'0'C^'0<0'0^>^^^>^•^^CD®®®®®0'0'0'0'0'0'0'0'0'0'0'0' 
^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  O  ^  ^  ^  O  4  ^  'O  *0  *0  *0  *0  'O  *0 


D  <  X  — 
t  ►-  O 
O  OC  ^ 
<  O  O  Of 
of  Z  »-  tu 


^(M<«^i/%a)(M'0'4r^(*^o®<04’  (*>^^1^000^  (>cDh><o  •^cD<Njmh-®0' 

®a)®®®0'0'00^(M<^J(*^•^u^<o^>®(^0'0^<^if*>4’u^'0'6r^^-®®(^(^0'(^0' 

^h“f*“^^^^®flOOO®OOCOCO®CO®CO®CO^O*^^^ff*^ff*^^^^^^^^^ 

^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  O  ^  ^  ^  O  ^  ^  ^  ^  ^  *0  ^  ^  O  *0  ^  *o  *o  *o  *o  *0 


«f'^<«1O'<9’®li^'t^•i^•<^0'f*^®®'^U^(^4.l/^O®^*’^0'<M(M0'U^4«^^>^(^J0'O®U^0D(M 

irNu^iA^oa)(>o<M'4’<o®^(*>'6®^'^r<*o(Niii>®'^(*>tAoDO(Ni(*>tA<or^®0'0'0' 

^  ^  ^  ^  ^  ^  ^  ^  ^  *0  ^  <0  ^  <0  <0  *6  <o  <0  <0  <0  <0  <0  *6  <o  ^  *6  <o  <0  <0  >o  <0  *6  <o  <0  <0  *6  <o 


ao«MO^<r^cDO<^®^^tA<o®rM»«f^mirNO.®9'h>ao<>rO'^(M»4r«jKh>«00'«rO'u>«4 

<Nf«><i’<.^«49'ir>iA«Du^h>'0'0rgo®>0(Nh>u>®*4  0''®«0oo>0Kir>«i4r»h>m®<6O 

<^J4’<^JU^<^<DO9'•^f*^OC^•««eD•^*^rMU^•-«•t^tf^00rr»l|^O9‘•OOl/^'>tOU^<tOlr>^tQ 

m9'®C^<^C^C^OU^f^'4'^>■U^U^O®O'0^>^J(N^»00'^'0'>t^C^•>^0«U^(^OO<^'0r<^ 

»OU^•>^^<NO’>0<^4®<N’>0<^«Mf<^4'U^^t't^JO®tf^-«^*moD^<>^*^JU^a'^J'OC^-^’>#’^• 

9'9'9'0'9'®®p^r*'O<rktr(-«rf<>rN*^O0'®'0iA'^(N*^9'®'0iAf<^«^O®'0irvr>>-4 


o«^<*^m>o<M«^9'0(ro®tf>mmrah>orN'0^c>>^oof<^r»(>j(>ja'oo(>joooio>^mh- 

^•®9'<Ntf^OU^fMO'^-IC<4'’<>^U^^O^>00O(^Jf<^l^^^0®®(^(^®^•U^rN®'4'®(^Jtf^^r*• 

a*c*<7*oo»^*^oj<N<*>’#<rv»oh-ooa*orvim^m»or»-oot^o^rvim't^'*rtf>i/>^^^»o 

^C*^^Cy'ff*^^^<^0‘0‘0‘0‘<7*0‘0‘CT*0‘0‘0'CT*0‘0‘CT‘0'0‘0‘0‘(7*(7*(7*0*<7*0'(7*C7* 


3  <  X  — 
M  X  »-  O 
O  ®  M 

<  o  o  oc 

OC  Z  »-  Ul 


(^J^^<Nu^m^f^<^Jfno'<^<^J®®o•«.ro^•'0'•r»^<^JC^u^o<^JCN®<^J<^J®c^^-(^^-o®o 

®®0'a*o-^-^<>jrvimm4'-4'tr\mm'0'0  '0'0>0 

'0'0'0'0'0'0'0>0'0'0'0'0'0'0'0<0'0'0'0h«h«h»h“h“r^h»f*“f*“r^h“f*“h“r^f*“h“h“r^ 
o*  o*  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^ 


0«^'00(rv«4®r»9'fn«4®®rvirg®o(^0(>J(>j®(r-<.ror*-h-'T'0<M'0'Oh-h--<>rr»'^ 

m  in  '0'0i0^»0'0p-h-h-h-r^h-h-h-h-h- 
^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^ 


omotnotnotrtotnotnou^oirtoi/^omotnomotnomomomomomo 

<^Jl^^^■o<^^tf'^•o<^J«f»^-o<N‘f'^■ocMlr^^-o<Mm^-o<M^f^^-o^^lm^-o<M^n^-o 

.^^^^<>i<>giM(>giv>mfnm^«r<Tsrmmmin^«0'0'0h-h-h-h-®®®®0' 


tA 

ec 

Ui 

Ui 

X 

o 


o 

< 

ac 


O 


z 

o 


</)  Ui  Ui 

3  X  ^ 

M  ^  ^ 

o  < 

<  Ik  ee 

X  o  < 


</) 
2  3 

Ui  O 

X  < 


^  Ui  < 
3  ^  X  — 
M  X  I-  o 
o  0^  M 

<0  0  ec 

oc  2  ^  Ui 


t/)  Ui  < 
3  X  « 
I-  O 
O  M 

<  u.  « 

0^  O  Ui 


o 

< 

ec. 


o 

(^ 

lA 

o 

h- 

o 

® 

O' 

O' 

iA 

iA 

O 

O' 

4- 

fA 

iA 

A* 

« 

O 

O 

4 

® 

4 

A- 

Ai 

4 

AI 

IA 

fA 

A- 

4 

IA  *4 

® 

h> 

® 

ru 

m 

O' 

® 

lA 

® 

o 

m 

o 

iA 

O 

O 

IA 

IA 

A* 

A- 

A- 

4 

AI 

O' 

® 

IA 

A» 

4 

4 

A» 

O 

IA 

4 

AI 

A-  O 

(A 

o 

fU 

fA 

h- 

Al 

® 

h- 

iA 

IA 

O' 

O' 

fA 

® 

fA 

Ai 

AI 

4 

o 

fA 

O 

A- 

Ai 

O' 

fA 

4 

O 

A» 

IA  O 

(A 

<M 

o 

m 

O' 

O 

A- 

Ai 

A- 

IA 

® 

® 

A- 

O 

o 

A* 

iA 

AI 

o 

O 

O' 

o 

® 

O 

O' 

AI 

iA 

O 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

9 

• 

• 

• 

9 

• 

• 

9 

9  9 

® 

(M 

<M 

® 

A4 

«A 

<# 

® 

o 

AI 

O* 

A- 

® 

4 

IA 

A- 

fA 

fA 

A- 

® 

fA 

4 

4 

w* 

IA 

A- 

O' 

r» 

h- 

in 

fu 

® 

Ai 

« 

® 

Ai 

® 

iA 

O' 

o 

o 

O' 

4 

fA 

O' 

IA 

O 

4 

® 

AI 

IA 

® 

fA 

IA 

A» 

m 

m 

m 

Ai 

Ai 

o 

c* 

® 

A- 

IA 

iA 

Ai 

o 

® 

A- 

IA 

iA 

® 

4 

4 

O' 

A» 

4 

w* 

O' 

4 

fA 

® 

IA 

ru 

lA 

lA 

lA 

«A 

IA 

IA 

IA 

IA 

4 

4 

4 

4 

4 

iA 

fA 

fA 

fA 

AI 

Ai 

AI 

AI 

•M 

O'  -t  >o 
^  'O 
^  «r  00 
00  <^  ^ 


<M  <^  »^  < 

o  o  « 
<^  u>  ^  < 
u>  ^  00  r 


^  ^  h- 

u>  u>  u> 

fO  fO  ^ 
^  'O 


r-  00  00  < 

U>  U>  lA  I 
m  r*>  f 
^  ^  « 


4  fA  4 

4 

Ai  4  O' 

® 

4  AI 

fA 

4 

O' 

•4 

® 

4 

4 

fA 

O  4 

A 

O' 

4  A  lA 

IA 

4  O' 

O' 

fA 

4 

A 

O 

® 

A 

4  4 

iA 

4  Ai  O 

A 

iA 

A  4  O 

A 

4  IA 

4 

AI 

A 

® 

® 

4 

O'  IA  Ai  O' 

A- 

IA  4  iA 

O 

®  4  4 

O 

4  ^ 

IA 

® 

O' 

O' 

® 

O  Ai 

fA  IA  A 

O' 

6  Ai 

fA 

4  IA 

4 

A 

® 

O' 

O' 

O'  O 

A- 

AAA 

®  ® 

®  ®  ® 

® 

O'  O' 

O' 

O' 

O' 

O' 

O' 

O' 

O' 

O' 

o  o 

fA 

iA 

fA  fA  fA 

fA 

fA  fA 

fA 

fA  fA 

fA 

^0  ^0  ^0  ^0 

4  4  4 

4 

4 

4  4  4 

4 

4  4 

4 

4 

4 

4 

4 

4 

4 

4 

4  4 

im«#<M^f*»f*«OO«#^>0<A00*^<^^'#>0U><AlA00lA^f*-*^00«0AI<Ah»^U><0 


o  r*->  00  ( 

ooooo«.,..  -_..  _ 

M»«»0000«^l*^l*^^OOf^^O<AOt^^^^-'----  -•  _-  —  _ 

mmu^f^<^<A^-^ioo<♦^<^^“U^^^ic^i^ooo^oor*•«o»r<^or^<^<^'^ooNu^^-oooo 

«coaoflOoo<^<^oo^<^i^ilA^u^^^«oo<^oo^^M<A«rll^*o*0^-^-ooooo^002 

f«i.f«»f«>f».fft.r**f^OOOOOOOOOOOOODODOOOOaOOD(^0*0'0*0'9‘0'09'9'9‘9*9*9*9‘9*9*9* 

mfnf*>miA<Am<A<AmiAmmrArAiAmrAmmrArArAm<A<ArAfAfnfr>mmfAtnfntntn 

0^'0^^^^'0^'0^^0^'000'0'0<4^'0^^'0«0'0'0<0'0<0«0'0'0'0'0'0 


^^f*»^^^lAf*»^OAif*\00OOCAO^f****''0Ai^OlAlArA^>0h“lA^^“*000^'>^ 

OOOf*>f^U^^MOO‘<^'0>fi'^<*'^^'^00^“0'OU^>#f*^OlA^^OOO•lAlAll^«OlA^•OAif'i'0 

tf>oo^-#®iA^mm<Mooo^«4‘mO'rrArA*4‘rA5KKr^^^®a*^m£^r^O'jO® 

f»>«AfAf*%<A<A®«AfAiAiAiAfA<A<A<ArArA<A<A<ArArAiArA<AmrAmrAiAmiAmmfom 


lA  Ui  Ui 

3  X  .J 
•-*  ^  «J 
O  < 
<  u.  ® 
a:  o  < 


\A 
Z  3 

<  M 

Ui  a 

X  < 


</)  U  Ui  < 
3  <  X  « 
^  X  I-  o 

a  ®  M 
<  o  a  ® 

®  2  *-  Ui 
X 


1/)  Ui  -< 
3  X  ^ 
M  ^  o 
O  M 

<  u.  ® 

®  O  Ui 

X 


IA  ^  A-  0^ 

4  A  IA  A  O  lA 

® 

O 

o 

O  4 

4  O 

pi4 

4  IA  O  IA  IA 

4 

IA 

o 

o 

IA 

IA 

pM 

A 

4 

® 

fA 

•M 

pri 

A  O 

4  IA 

pM  o  fA  AI  4  AI 

4 

® 

i^ 

®  A 

A  IA 

o 

fA  A  lA  O  4 

PA 

® 

o 

o 

AI 

P^ 

o 

® 

IA 

fA 

iA 

O 

O  4 

Ai  Ai 

O  fM  lA  A4 

fA 

p4 

4  4 

P^  o 

IA 

O  A  O  i'i  O 

O 

4 

prf 

pM 

4 

A 

IA 

A 

A 

iA 

» 

4 

fA  IA 

fA  4 

4  O'  O  A  4 

O 

®  «-« 

AI  Psi 

IA 

AI  4  4  »  4 

IA 

4 

o 

4 

O 

IA 

4 

O 

IA 

4 

O 

IA 

4 

•  • 

•  ••••• 

• 

• 

• 

•  • 

•  9 

• 

•  •  •  •  • 

• 

• 

9 

9 

• 

• 

• 

• 

• 

• 

9 

• 

• 

iA  O' 

®  <^ 

i<l  O  CA  p^  lA  fA 

4 

® 

IA 

IA  O 

®  o 

4 

A  AI  PsI  A  ® 

4 

4 

4 

o 

pi4 

O* 

IA 

o 

O 

O 

o 

4 

fA 

4  lA 

4  AI 

O  4  AI  ®  AI  4 

O 

Pi4 

fA 

4  IA 

4  4 

AI 

O  ®  IA  p4  A 

iA 

® 

iA 

A 

PSl 

IA 

O 

AI 

4 

o 

p4 

4 

A 

O'  o 

O  O 

O  ®  ®  A  A  4 

IA 

IA 

4 

iA  AI 

o 

O 

®  4  IA  4  AI 

■•4 

O 

® 

4 

lA 

fA 

o 

® 

4 

in 

iA 

iA  iA 

fA  iA 

fA  fA  ^  fA  ^  fA 

fA 

fA 

fA 

fA  fA 

fA  iA 

AI 

AI  AI  Ai  AI  Pi 

AI 

p4 

p4 

Q'Q'O'O'O'O'O'O'O'O'O'f^OO'O'O'OOOO'OOO'OOO'OO'O'OOO'OOO'O'O' 

jAf<ifnmfnfnfn«AfnfniA»AfomfnfnfAf<>fOfnfOfrifOfnfOfnfnfnfnfnfAfAfnfOfnc<>fA 


tf>^h-tf>mp^fvr^®O^Of>i-**®®'»’0't^^^-^r^po^i^^w^cjco^®flO»Aj^ 

fnfn-^w>f»»OfMiA®rg«o-M®-M®f>i®«t^«>-^r-rgw»A®f>ir^^«rr^OfM'tu>.0® 

Q>Q>0>Q'9>Q'0‘0‘^f^OO‘00'000'0'0'000000000'0'0'0'0'0'000‘0‘ 

fnmfnfnfnfnfnfnmmmfAfo»AfnfAfAfOfr><^<^<^<AfAfnmfnfnmfntnfnfn<AfntOfn 


^Q'Q'Q'0'0'^0'^^00'0>0000000000000000000'00'00  0 
fOfA<n<Afn«AfOfnfAfOfOfOfOfOfAfAfA<A<AfA«AfAfAfAf<><^<AtotoMtoM<AtoM<Af«> 


OU>OU>OIAOlAOIAOU>OU>OU>OU>OlAOU>OlAOU>OU>OU>OIAO(AOU>0 


7 


c)  Record  the  observations  in  latitude  and  longitude,  apply  a 
transformation  to  obtain  rectangular  coordinates,  and  then 
perform  the  claculations  using  the  plane  formulae  for  the 
average  and  variance. 

Procedure  (a)  has  the  disadvantage  of  being  more  complicated.  A 
sufficiently  small  portion  of  the  earth's  surface  can  be  considered  a 
plane  and  the  additional  complication  introduced  by  the  use  of  spherical 
versions  of  the  statistical  formulae  may  not  be  warrented.  Somewhat 
similar  problems  have  been  investigated  in  the  field  of  land  surveying 
and  are  reported  in  most  works  on  geodesy.  Procedures  (b)  and  (c), 
above,  are  mathematically  equivalent  since  maps  are  made  by  transforming 
latitude  and  longitude  to  plane  coordinates  via  a  map  projection.  Hence 
a  study  of  the  numerical  differences  between  computations  on  a  plane  and 
on  the  earth  becomes  a  study  of  map  projection  distortions. 

The  official  map  producing  agencies  of  the  various  countries  of  the 
world  have  recognized  the  advantages  of  rectangular  coordinates  for  local 
purposes  and  save  the  map  user  the  trouble  of  assigning  his  own  system 
of  rectangular  coordinates.  They  do  this  by  publishing  maps  which  have 
the  official  plane  coordinates  printed  directly  on  the  maps.  Two  map 
projection  systems  of  this  type  are  available  in  the  United  States,  and 
comparable  systems  exist  in  most  other  countries  of  the  world.  The  use 
of  these  systems  is  not  restricted  to  calculations;  they  might  also  be 
used  to  record  and  index  information  in  terms  of  the  plane  coordinates, 
perhaps  scaled  from  topographic  maps.  The  systems  now  available  have 
several  features  in  common.  The  coordinates  are  usually  given  as  rectan¬ 
gular  coordinates,  often  chosen  so  that  all  values  are  positive.  More 
importantly,  the  errors  in  computing  as  though  the  earth  were  a  plane 
disk  can  be  evaluated.  This  implies  that  the  region  within  which  one 
can  perform  plane  computations  with  a  specified  level  of  precision  can 
be  defined  on  an  a  priori  basis.  If  the  allowable  error  is  small,  the 
region  must  be  small  or  several  map  projection  systems  (called  zones) 
must  be  used  within  the  region.  In  the  latter  event  conversion  between 
zones  may  be  required.  This  conversion  may  be  directly  from  zone  to 
zone  or  may  involve  reconversion  to  geodetic  coordinates  as  an  intermediate 
step.  There  are  certain  advantages  in  using  a  conformal  map  projection 
for  such  a  system  since  the  scale  errors  are  then  independant  of  direction 
and  a  scale  factor  can  be  applied  to  improve  the  accuracy  of  short  lengths  . 
The  two  systems  employed  in  the  United  States  are: 

1)  The  State  Plane  Coordinate  System:  This  system  comprises 
approximately  120  zones  covering  the  entire  United  States,  with  the 
orientation  toward  individual  states.  The  accuracy  within  each  zone  is 
one  part  in  10,000.  The  larger  states  therefore  require  several  zones. 

The  zones  overlap, with  boundaries  between  zones  lying  along  minor  civil 
divisions  (usually  counties).  The  Lambert  Conformal  Conical  projection 
and  the  Transverse  Mercator  projection  are  employed  (with  only  one  exception) , 
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depending  on  the  shape  of  the  individual  states.  This  system  is  adtnira  y 
suited  to  the  needs  of  the  local  land  surveyor  and  has  been  officially 
adopted  by  many  local  governmental  units.  In  many  states  it  has  legal 
status,  is  used  for  land  ownership,  and  appears  on  large  scale  maps.  Con¬ 
version  tables  are  available  and  simple  to  use  for  any  particular  zone. 
Conversion  between  zones,  and  especially  between  states,  is  somewhat  more 
inconvenient.  The  location  of  the  zones  occasionally  is  awkward.  In 
Washington  state,  for  example,  the  two  merging  metropolitan  areas  of 
Seattle  and  Tacoma  each  lie  in  a  separate  zone.  The  system  of  State  Plan 
Coordinates  appears  on  all  recent  U.S.  geological  Survey  topographic  maps. 

2)  The  Transverse  Mercator  System: 

Known  as  the  Universal  Transverse  Mercator  grid  system  (UTM) 
this  system  is  employed  by  the  U.S.  Arny.  The  UTM  grid  extends 
to  eighty  degrees  north  and  south  latitude,  beyond  which  a 
Polar  Stereographic  grid  is  employed.  The  UTM  grid  extends  around 
the  world  in  sixty  north-south  zones,  each  covering  six  degrees 
of  longitude  with  an  overlap  of  one  half  degree.  The  accuracy 
within  each  zone  is  one  part  in  2500.  Since  different  areas  of 
the  world  are  based  on  distinct  ellipsoidal  datums,  separate 
tables  are  required  for  various  parts  of  the  world.  Procedures 
are  available  for  converting  directly  from  one  zone  to  adjacent 
zones.  The  zonal  nature  of  the  system  is  occasionally  inconvenient. 

An  alphanumeric  partitioning  of  areas  is  available  in  the  system. 

The  UTM  grid  appears  on  all  Army  Map  Service  topographic  maps, 
some  foreign  maps,  and  on  all  recent  U.S.  Geological  Survey 
topographic  maps. 

Both  of  the  foregoing  systems  have  several  advantages.  They  '^^n 
be  employed  for  virtually  all  computations  without  serious  error.  Further, 
any  information  recorded  in  either  of  these  systems  can  be  related  ..o 
genetic  coordinates  and  hence  to  information  collected  anywhere  else  in 
fhe  world.  Also,  these  coordinates  are  already  shown  on 
and  most  photogrammetric  firms  are  sufficiently  familiar  with  these 
systems  to  add  them  to  aerial  photographic  or  maps  compiled  by 
trie  methods.  The  disadvantages  of  these  systems  stem  largely  from  their 

can  be  circumvented  in  several  ways. 

When  a  map  projection  system  is  to  be  used  soley  for  computational 
purposes,  and  not  necessarily  to  be  indicated  c.i  published  "«P®>  ^ 
choice  of  a  particular  projection  depends  on  the  type  of  computation 
contemplated.  The  systems  cited  above  are  so  ^ 

stated  level  of  accuracy  for  virtually  all  computations.  This  is 
a  restriction  which  narrows  the  range  of  suitable  projections  and  resul 
in  projection  which  require  fairly  involved  computations.  For  a  gi 
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problem  there  may  be  a  specific  projection  which  is  computationally 
much  simpler  but  which  yields  results  which  are  of  equal  accuracy.  For 
example,  a  problem  which  requires  interpolation  between  two  points 
on  a  sphere  might  be  attacked  by  using  the  gnomonic  projection  (see 
Appendix)  since  all  great  circles  are  straight  lines  on  this  projection; 
linear  interpolation  in  gnomonic  coordinates  will  yield  a  point  lying 
on  the  arc  connecting  the  two  given  points.  Similarly,  problems  involving 
circles  on  a  sphere  may  be  attacked  using  the  stereographic  projection. 

In  other  situations  computational  simplicity  and  speed  may  be  more 
important  than  a  few  tens  of  meters  of  accuracy,  Kao,  for  example,  has 
recently  shown  that  the  geometric  (perspective)  projections  are  especially 
well  suited  for  calculation  by  digital  computer,  particularly  when  large 
amounts  of  locational  information  are  required  within  fractions  of 
a  second  (i.e,,  in  real  time  problems).  Clearly  the  chaice  depends  on 
the  nature  of  the  problems  and  the  volume  of  the  information  to  be  processed. 
Computer  calclation  of  distance  and  direction  on  a  sphere  (or  ellipsoid) 
may,  in  many  instances,  be  easier  than  attempting  to  convert  to  plane  co¬ 
ordinates,  On  the  other  hand  a  more  complicated  problem,  as  for  example, 
occurs  in  weather  prediction,  may  advantageously  be  solved  by  the  use  of 
an  appropriate  map  projection.  In  this  instance  the  problem  is  to  con¬ 
struct  contour-type  maps  of  the,:  entire  northern  hemisphere  from  informa 
tion  received  from  locations  scattered  within  this  region.  Rather  than 
attempting  to  solve  the  contour  interpolation  problem  on  a  sphere,  the 
Weather  Bureau  employs  stereographic  map  projection  coordinates  with 
a  local  correction  for  the  projection  distortion  and  solves  the  problem 
in  plane  coordinates. 

If  one  has  information  recorded  in  latitude  and  longitude  simple 
conversions  to  map  projection  coordiantes  are  available.  For  example, 
one  can  pretend  that  these  are  already  the  ordinate  and  abscissa 
of  a  plane  coordinate  system.  The  resulting  projection  is  known  as  the 
square  projection.  Computations  performed  in  this  manner  will  differ 
from  the  true  values  by  amounts  which  depend  on  the  size  of  the  region 
and  on  the  latitude.  Another  simple,  but  slightly  better,  conversion  is 
to  multiply  all  the  abscissas  (longitudes)  by  a  constant  equal  to  the 
cosine  of  the  average  latitude  of  the  region  in  question  (the  square 
projection  with  a  standard  parallel;  also  known  as  the  rectangular  pro¬ 
jection),  Such  a  procedure  might  for  example,  by  employed  in  urban 
analysis,  depending  on  the  size  of  the  area.  Another  alternative  would 
be  to  transform  to  rectangular  coordinates  by  converting  all  values 
into  distances  north  and  east  (that  is,  measured  along  a  parallel) 
from  some  arbitrary  point  within  the  region.  This  yields  the  sinu¬ 
soidal  projection.  The  equations  for  all  of  the  above  projections  are 
extremely  simple.  Somewhat  more  refined,  but  also  more  complicated,  sol¬ 
utions  take  into  cosideration  the  shape  of  the  area  of  concern,^  Albers 
equal  area  conical  projection  with  standard  parallels  at  29  30  N  and  45 

30'N,  and  Lambert's  conformal  conical  projection  with  standard  parallels 
at  33°N  and  45°N,  for  example,  are  two  systems  which  might  be  suitable 
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for  the  continental  United  States.  The  distance  error  in  computing 
with  these  latter  systems  is  not  likely  to  exceed  fifty  miles. 

The  use  of  latitude  and  longitude,  while  advantagdl(f  from  the  point 
of  view  of  long  run  national  needs,  entails  some  local  difficulties.  In 
the  process  of  recording  it  may  be  necessary  to  interpolate  between  curved 
lines,  and  the  system  of  minutes  and  seconds  is  awkward  (Decimal 
degrees  are  more  convenient).  The  complete  number  of  digits  required 
to  specify  a  given  location  in  its  world  context  is  excessively  large 
for  local  use,  and  the  north-south  and  east-west  designation  is  often 
superflouous  (a  mathematical  convenience  is  obtained  if  south  latitudes 
and  west  longitudes  are  considered  negative).  Finally,  and  perhaps 
most  important,  it  is  often  difficult  to  determine  the  latitude  and 
longitude  of  a  particular  spot. 

The  direct  recording  and  storage  of  geographical  information  in 
terms  of  rectangular  coordinates  circumvents  some  of  these  difficulties, 
but  introduces  others.  The  majority  of  the  electro  -  mechanical  data 
reduction  devices  (specifically,  coordinate  readers)  which  are  now  on 
the  market  utilize  rectangular  coordinates.  These  instruments  reduce 
the  teduim  of  coordinate  reading,  even  when  the  desired  result  is 
latitude  and  longitude  coordinates.  In  this  case  the  inverse  map  pro¬ 
jection  equations  are  required.  Curiously,  these  are  not  widely  avail¬ 
able  in  the  literature  on  the  subject  of  map  projections  (with  a  few 
exceptions)  since  the  previous  technology  prohibited  their  extensive 
use. 


If  the  objective  of  the  study  does  not  include  subsequent  conversion 
to  latitude  and  longitude,  a  convenient  procedure  is  to  draw  arbitrary 
rectangular  (or  polar)  coordinates  on  whatever  maps  or  aerial  photographs 
are  available.  One  advantage  is  that  this  can  be  done  by  persons  with 
no  training  and  with  virtually  no  intellectual  effort  or  financial 
expenditure.  When  the  map  used  is  accurate  and  at  a  "sufficiently  large" 
scale  these  arbitrary  coordinates  may  be  employed  as  are  the  map  projection 
coordinates  discussed  above.  If  the  information  collected  has  no  permanent 
value,  this  procedure  is  perfectly  satisfactory. 

A  disadvantage  is  that  the  errors  introduced  are  not  known.  The  limits 
within  which  a  certain  level  of  accuracy  obtains  is  uncertain  an  one  never 
knows  whether  the  system  can  be  extended  to  include  a  neighboring 
territory.  A  second  major  disadvantage  is  that  it  may  not  be  possible 
to  use  information  collected  for  oie  study  in  a  second  study  which  either 
(a)  encompasses  a  larger  area  than  the  original  study,  or  (b)  which  is 
subsequent  in  time  to  the  original  study,  especially  if  the  original  map 
has  been  lost,  or  (c)  which  requires  a  higher  level  of  accuracy  than  the 
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the  original  study.  One  can  imagine  the  difficulty  of  analyzing  the 
greater  metropolitan  area  of  Kansas  City  if  Kansas  City,  Missouri  and 
Kansas  City,  Kansas,  used  two  different  and  unrelated  grid  systems. 

Or  if  each  bureau  of  a  city  government  employed  a  distinct  system  of 
coordinates.  The  actual  occurance  of  situations  of  this  very  nature 
in  the  field  of  civil  engineering  is  what  gave  impetus  to  the  estab¬ 
lishment  of  the  system  of  State  Plane  Coordinates  by  the  U.S.  Coast  and 
Geodetic  Survey  in  the  1930 's. 

Conversion  between  arbitrary  map  coordinates  can  be  effected  with 
relative  ease  if  the  relation  between  the  two  systems  is  known,  or  if 
both  systems  are  related  to  latitude  and  longitude  by  known  inverse 
equations.  If  the  relation  betv?een  systems  is  not  known  it  is  theoreti¬ 
cally  possible  to  estimate  the  relation  if  the  coordinates  of  a  suffi¬ 
cient  number  of  points  are  accurately  known  in  both  systems  (see 
Appendix).  Such  conversions  may  occasionally  be  required  but  are 
expensive. 

A  final  distinction  should  be  made  between  coordinates  and  areas. 
Coordinates  describe  points,  not  areas,  and  one  must  distinguish  be¬ 
tween  an  real  recording  unit  such  as  a  census  tract  and  between  the 
coordinate  system  used  to  pinpoint  some  centroid  taken  to  represent 
that  areal  unit.  Areal  information  recording  units  are  extremely 
numerous  and  differ  widely  in  size  and  shape.  As  a  consequence  it 
is  often  necessary  to  convert  from  one  areal  unit  (e.g.  census  tract) 
to  other  areal  units  (school  district,  political  precint,  and  so  on). 
These  areal  conversions  differ  somewhat  from  the  coordinate  conversions 
discussed  in  this  report.  In  general,  specification  of  the  areal  bound¬ 
aries  must  be  included  in  the  mathematical  conversion  statements. 

There  are  then  again  several  procedures,  of  varying  accuracy  and 
complexity,  which  may  be  employed  for  the  conversions. 
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APPENDIX  I 


CONVERSION  FROM  THE  PUBLIC  LAND  SURVEY 
SYSTEM  TO  LATITUDE  AND  LONGITUDE 


The  simplest  conversion  begins  with  a  procedure  which  assumes  that 
the  Public  Land  Survey  conforms  to  the  exact  specifications  upon  which 
it  is  based.  The  system,  as  is  well  known,  does  not  conform  to  these 
specifications,  for  a  number  of  reasons  including  measurement  errors 
unavoidable  in  any  empirical  work  and  a  certain  laxity  of  supervision 
during  the  establishment  of  the  system.  For  conversion  into  latitude 
and  longitude  the  following  notation  is  convenient: 


>i 


a 


e 


Mi 


T 

R 

^n 

Se 

On 

Qe 

0 


is  an  index  to  indicate  the  initial  point  of  the  survey.  It  is 
necessary  to  distinguish  at  least  37  initial  points  in  the  Western 
United  States. 

is  the  latitude  of  the  i^^  base  parallel. 

is  the  longitude  of  the  i^^  base  meridian,  with  west  longitudes 
negative. 

is  the  equatorial  radius  of  the  ellipsoid  taken  to  represent  the 
earth.  For  the  Clarke  Ellipsoid  of  1866,  a= 3963. 2257  miles. 

is  the  eccentricity  of  the  ellipsoid  taken  to  represent  the  earth. 
For  the  Clarke  Ellipsoid  of  1866  e  =  0.0822718542. 

is  the  radius  of  the  meridian  at  the  i*^^  initial  point.  Mj^  is  given 

by  , 

a(l  -  e^) 

<l  -  .2  sln^ 

is  the  township  number  of  the  location  in  question,  with  north 
townships  taken  as  positive  and  south  townships  taken  as  negative. 

is  the  range  number  of  the  location  in  question,  with  east  ranges 
taken  as  positive  and  west  ranges  taken  as  negative. 

is  the  northing  of  the  section  in  question,  with  the  sign  convention 
as  above. 

is  the  section  easting,  with  the  sign  convention  as  above, 
is  the  quarter  section  northing,  with  signs  as  above, 
is  the  quarter  section  easting. 

is  the  latitude  (to  be  found)  of  the  location  in  question. 


is  the  radius  of  curvature  perpendicular  to  the  meridian  at 
latitude  0; 


The  necessary  equations  are  then: 

6  T  T  3  +  ^ 


6  R  T  3  +  Sg  -K  Qg 
cos  0 

The  formulae  are  established  by  observing  that  the  center  of  the 
township  in  question  should  be  six  miles  times  the  number  of  the  township 
north  (south)  of  the  base  parallel,  minus  three  miles  to  obtain  the 
center  of  the  township.  The  section  northing  and  easting  give  the 
distance  of  the  center  of  the  section  from  the  center  of  the  township, 
and  the  quarter-section  northing  and  easting  give  the  distance  of  the 
center  of  the  quarter-section  from  the  center  of  the  section.  For  the 
SW  1/4,  Sec.  25,  T  5  N,  R  17  E,.  one  should  have  for  example,  that  the 
center  of  the  township  is  27  miles  (5x6-3)  north  of  the  initial  point. 

S  is  -1.5  miles  and  Q  is  -0.25  miles.  The  total  distance  in  the 
north-south  direction  from  the  initial  point  should  therefore  be  +  25.25 
miles.  This  distance  must  then  be  converted  to  the  appropriate  number  of 
degrees  and  added  to  the  latitude  of  the  initial  point.  A  further  refine¬ 
ment,  though  hardly  necessary,  would  be  to  iterate  on  the  latitude  obtained 
in  the  first  step  in  order  to  adjust  the  meridional  radius  employed  in  the 
computation.  Determination  of  the  longitude  is  similar  but  slightly 
more  difficult  since  the  distances  are  measured  along  a  parallel  (a 
loxodrome,  not  a  great  circle  or  geodesic),  whose  radius  varies  with  the 
latitude. 

In  programming  the  outlined  procedure  for  a  digital  computer  it  is 
simplest  to  employ  radians  instead  of  degrees  and  to  store  a  table  of 
Sn,  Sg,  0^,  Ofj,  0i.  and  Ai.  The  computer  can  perform  the  assignment 
to  the  correct  initial  point  by  letter  for  letter  examination  of  the  name 
of  the  principal  meridian.  A  convention  is  necessary  to  distinguish 
the  two  different  initial  points  employed  for  the  Fourth  Principal  Meridian. 
The  method  detailed  assigns  latitude  and  longitude  (to  about  the  nearest 
1/4  mile)  on  the  assumption  that  the  Public  Land  Survey  designations  are 
where  they  should  be.  Of  course  they  are  not  exactly  there:  the  legal 
strategy  is  to  assign  to  the  actual  locations  a  status  of  incontestable 
correctness,  irrespective  of  any  errors  which  may  have  been  introduced 


0  '  0i  ^ 

and 


during  the  survey.  To  adjust  the  calculated  values  to  conform  to  their 
legal  positions  requires  detailed  historical  and  empirical  corrections, 
and  can  be  quite  tedious.  For  n»any  research  purposes,  hov;ever,  such  a 
refinement  may  not  be  necessary.  To  obtain  an  order-of-magnitude  estimate 
of  the  discrepancies,  the  actual  latitude  and  longitude  (as  recorded  on 
large  scale  topographic  maps)  of  a  scattered  set  of  locations  have  been 
compared  with  the  computed  values.  For  a  selection  of  74  points  within 
the  State  of  Michigan  the  errors  are  as  follows: 

Distribution  of  Errors:  (N  =74) 

Mean:  2.849  miles 

Standard  deviation:  1.827  miles 


Maximum:  9.339  miles 

6074  of  the  errors  are  less  than  2  miles 


937o  of  the  errors  are  less  than  5.5  miles 


The  directional  errors  appear  evenly  distributed  in  all  directions. 
A  random  selection  of  points  (N  -  25)  from  other  states  indicates  that 
the  errors  are  quite  comparable  and  of  the  same  order  of  magnitude.  A 
sample  computation  is  as  follows: 


observed  location: 
calculated  Lat/Lon: 
observed  Lat/Lon: 
difference : 


Svl  1/4,  Sec.  28,  T  2  S,  R  6  H,  Michigan  Meridian 


o 

CM 

16' 

07' 

N, 

CO 

o 

44' 

08"  i; 

42° 

17' 

10'' 

N, 

83° 

44' 

49"  W 

1' 

'03" 

41'' 

difference  in  miles:  2.  5G 

direction  of  difference;  154.34°  (E  of  N) 

The  method  given  above  does  not  include  an  adjustment  for  the  con¬ 
vergence  of  the  meridians.  Since  the  edges  of  the  ranges  run  due  north 
and  south,  the  ranges  become  narrower  as  the  meridians  converge.  To  adjust 
for  this,  standard  parallels  are  established  every  twenty-four  miles  north 
and  south  of  the  base  parallel.  The  ranges  are  again  made  six  miles  wide  at 
these  standard  parallels.  The  system  thus  is  self  correcting  every  twenty- 
four  miles.  The  order  of  magnitude  of  the  difference  in  width  of  ranges, 
separated  by  twenty-four  miles  in  a  north-south  direction,  can  be  established 
as  follows:  The  radius  of  the  parallel  at  45°N  latitude  is  2807,178  miles. 

At  a  distance  of  24  miles  north  of  45°N  it  is  approximately  2789.834  miles. 

The  east-west  width  of  the  northern  edge  of  the  range  24  miles  north  of 
the  45th  parallel  is  therefore  not  six  miles  but  0.037104  miles  (195.9  feet) 
less  than  six  miles.  On  this  basis  the  error  at  R  30  E,  an  extreme  value, 

would  be  1.86  miles.  Another  slight  error  is  introduced  by  the  topographic 

elevation,  since  the  radii  employed  apply  to  a  mean  sea  level  ellipsoid. 


Empirical  corrections  for  Michigan  vjould  need  to  Include  the  fact  that  the 
standard  parallels  are  60  (not  24)  miles  apart  (In  accord  with  the  surveying 
Instructions  In  force  at  the  time),  and  that  R  1  E  is  consistently  too 
narrow  from  T  1  N  to  T  20  N.  An  adjustment  for  these,  and  other,  system¬ 
atic  departures  could  be  Incorporated  Into  the  computer  program.  Conver¬ 
sion  of  the  Section,  Township,  and  Range  Information  to  latitude  and 
longitude  can  be  followed  by  conversion  to  map  projection  coordinates 
for  map  plotting  or  computational  purposes.  Direct  conversion  to 
Cartesian  map  coordinates  also  Is  possible  but  is  less  convenient  for 
the  entire  Western  United  States.  This  is  more  appropriate  for  operations 
restricted  to  a  limited  area,  e.g.,  one  individual  state. 
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APPENDIX  II 


MAP  PROJECTION  EQUATIONS 

The  following  list  gives  the  mathematical  rules  for  the  most  common 
map  projections  of  a  sphere.  The  following  notation  is  standard. 

(f  Latitude  of  a  point  whose  projection  coordinates  are  desired. 

^  Longitude  of  a  point  whose  projection  coordinates  are  desired. 

Abscissa  of  a  plane  cartesian  coordinate  system. 

'y  Ordinate  of  a  plane  cartesian  coordinate  system. 

r  Radial  distance  of  a  plane  polar  coordinate  system. 


Angular  direction  of  a  plane  polar  coordinate  system. 

Latitude  of  the  center  of  the  map;  either  the  point  of  "tangency", 
or  a  single  standard  parallel. 

Southerly  standard  parallel  for  projections  having  two  standard 
parallels. 


Northerly  standard  parallel  for  projections  having  two  standard 
parallels. 


^  Longitude  of  the  center  of  the  map;  either  the  point  of  "tangency", 
or  the  central  meridian. 

C.  The  constant  of  the  cone  for  conic  projections. 

The  radial  distance  from  the  origin  to  the  image  of  the  southerly 
*  standard  parallel  in  plane  polar  coordinates. 


Air  equations  are  giver  for  a  sphere  bf  unit  radius  (R  =  1)  and  all 
values  are  assumed  to  be  in  radians.  Conversion  to  scale  can  be  achieved 
by  multiplying  all  distances  by  the  appropriate  scale  factor.  North 
latitudes  and  east  longitudes  are  taken  to  be  positive,  i.e. 

-  jT 

-ir  f  /  ^ 

The  equations  are  given  in  their  most  commonly  applied  form.  The 
conical  projections,  for  example,  are  not  given  in  their  oblique  cases. 
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i 
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(1)  Albers'  equal  area  conic  projection  with  two  standard  parallels; 


C  -  ^  -A  ^//7 


r  = 


i/a 


c 


This  puts  the  origin  of  the  coordinates  somewhat  beyond  the  north 
pole,  which  is  rather  inconvenient.  The  origin  can  be  shifted  to  the 
intersection  of  the  southern  standard  parallel  with  the  central  meridian 
by  using 

X  =  ^  ^ 

y  =  V,  -  ^  ^ 


(2)  Azimuthal  equidistant  projection  *. 

r=  drc.  f  ^  Cos(x  -  J 

CoS  tp  Sin  CA-Xa")l 
Sm  >r 

The  origin  of  the  coordinates  is  at  ^  X  o 

(3)  Bonne's  Equal  Area  projection  t 

—  Lp  4-  ”fc<Xv'  "x 

^  -  Xo)  ^  ^  ^ 


9  =  Cure  Siv\ 


1  “* 

— 
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To  place  the  origin  at  the  intersection  of  the  standard  parallel  and 
the  central  meridian  use: 

r,  =  -(-(Xn  (  ^  "  Vo  ) 

X=  r  B 

Y  =  r,  -  r  ^  . 


(4)  Cassini  Projection 

X  -  Cl<C  ^  i  \n 


A'Vyx 


Kj  -  —  +-  (k.'TC. 


-x.)] 

'to.n  ^ 

C.O-S  Cx  ~  X  c  ^  J 


(5)  Gnomonlc  Projection: 


X  - 


C£>^  ^  (\  ~\  o') 


Y 


Sitn  vf  SiK  +  Costf  Cos  C.o^(X-  Xo) 

Sm  4?  Os  -  Sin  (Loh  ^  <L<n£  (\  - 

'f  If,,  ^  <LoS  f  CoS  (po  CoS  (X  -  Xo') 

(6)  Lambert's  azimuthal  equal  area  projection: 

t  ^  drc  If  •+  CexS  \p  CoS  Cos  (>.-  Xo)] 

r=  a 

tf  (^X~  X  cJ  ^ 


0  =  A.r<:  Sin 


Sm  C 


(7)  Landsert's  cylindrical  equal  area  projection: 

X  =•  X  -  Xo 
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Or,  with  a  standard  parallel; 

X  =  (  X  -  y  o  ')  Co  S  % 

\y  =  ^ 

(8)  Lambert's  conformal  conic  with  two  standard  parallels; 

Ift  Cog  —  Lh  (fz. _ 

Lvx  4-d.h  ( '  -  ~  ^  ^  ) 

Q-.  d  (\  -  X=) 

X  ^  r  ^lyy  9 

Y  =  r,  -  'T  Cos,  G 

(9)  Mercator's  conformal  cylindrical  projection; 


X  -  X  -  Xo 


(10)  Miller's  Cylindrical  projection: 

X  -  X  - 

Y_  /'«  ~h&H  [Tf  *'  ^ 
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(11)  Mollweide's  equal  area  elliptical  projection: 

Define  chi  by  ^  2.  ^  Siry  tp  ^ 

( A  -  >o)  (:!os  42 

Y  ^  V 

(12)  Orthographic  projection: 

X  =  SiA  4>  Cjos  4)„  -  (L<9S  f  ("a  -  Ac  ) 

\l  ^  (Los  f  (x- Xo  *) 

(13)  Polyconic  projection  (American  poly conic) : 

r-  f 

0  =  (X-  X.')  V" 

)<  ,  \r  I  h  0- 

■y  _  <P  —  (Lo  S  © 

Which  puts  the  origin  at  the  equator. 

(14)  Sinusoidal  equal  area  projection: 

y-  (x-Xo)  c-o-^  f 

y  =  ^ 

(15)  Square  projection: 

\  -  X  -  X= 

7  -  r 


or,  with  a  standard  parallel  (also  known  as  the  rectangular  projection); 

x*  (x  -X.') 

y-  ^ 


(16)  Stereographic  projection: 


X  = 


Co-"^  ^ 


^4-  ■*"  ‘f’  C^~  ^o  ) 


y  - 


vf  CiPSip,  -  dose/)  (Vq^  (>V-Xo) 

1  4-  ?.Ui  *^0  ^  Cps(  A  ~  Xo) 


(17)  Transverse  Mercator  projection. 

1  4  (!oS.  4^ 

4.  ~  C.<o^  ^  (.  ^  ~  Xo') 


y  --  ^  Lh 


'y  =  "^Ay*  |*4ii-vxy'  ^X  —  "X.*^ 


APPENDIX  III 


LEAST  SQUARES  CONVERSICW  FROM  ONE  SYSTEM  OF 
RECTANGULAR  COORDINATES  TO  ANOTHER 

Given  two  sets  of  coordinates  on  the  same  map,  with  a  minimum 
of  five  points  identified  in  both  systems  of  coordinates,  it  is  possible 
to  convert  the  coordinates  of  one  set  to  the  other  by  a  two-dimensional 
version  of  a  least-squares “line"  The  procedure  is  most  easily  effected 
using  complex  numbers. 

Let  x,y  be  one  set  of  coordinates  and  u,v  be  the  other  set,  and 
let  Wj  :  X  iy  and  Zj  r  u-e  iv,  where  i^  =  -1,  be  the  complex  numbers 
representing  the  i^^^  point.  The  objective  is  then  to  find  the  complex 
constants  A  =  aj^4-  ia2  and  B  =  bj^+  ib2  in  the  equation  W  :  A  +•  BZ  such 
that  the  squared  residual 

2  1  Wj  -  w- 1 

j'-' 

is  a  minimum.  The  normal  equations  are  readily  obtained  by  differentia¬ 
tion.  The  equation  can  be  rewritten  as  a  pair  of  transformation  equations 
by  separating  the  real  and  imaginary  parts,  viz: 

A  A 

Re  (W)  =  X  =  3]^  +•  b]^  u  ^  b2v  • 

A  A 

Im  (W)  r  =  32  ♦-  b2  u  -  b^v' 

where  x  and  y  are  the  estimates  of  the  x,y  coordinates  as  obtained  from 
the  known  u,v  coordinates.  The  standard  error,  etc.,  of  the  estimate 
can  be  obtained  in  a  manner  analogous  to  that  employed  for  ordinary 
least  squares  procedures. 


A  similar,  but  considerably  more  complicated,  procedure  must  be 
employed  if  the  two  sets  of  coordinates  do  not  come  from  the  same  map, 
or  if  the  relation  to  latitude  and  longitude  is  to  be  estimated,  or 
if  an  attempt  is  made  to  determine  the  map  projection  of  an  arbitrary 


map. 
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R  RETURNS  DISTANCE  IN  KILOMETERS 
R  RETURNS  DIRECTION  IN  DECIMAL  DEGREES 
R  ACIC  TR  80.  PAGES  41-47. 

R  NECESSARY  CONSTANTS 
VECTOR  VALUES  P I =3 1 4 1 5 9265 E-8 
VECTOR  VALUES  T P I =6283 1 8 5 3 1 E-8 
VECTOR  VALUES  AR AD=63 782064E-4 
VECTOR  VALUES  BR AD =6 3 565 838E-4 
VECTOR  values  BO VRA=9 966099 247 E-1 0 
VECTOR  VALUES  VK 1 = 2 3 562 1 8428E-7 
VECTOR  VALUES  VK2=6956258069E-5 
VECTOR  VALUES  VK 3 =4986428206 E-9 
VECTOR  VALUES  VK 4= -40 10886986E -1 0 
VECTOR  VALUES  VK 5=-79945 5 6 507E-1 0 
VECTOR  VALUES  VK6=3986428206E-9 
VECTOR  VALUES  E1=17036962E-10 
VECTOR  VALUES  E2=21769E-10 
VECTOR  values  E3=29026E-10 
VECTOR  VALUES  E4=3628E-10 
VECTOR  values  RAD=174532925E-10 
R  BEGIN  COMPUTATION 
ENTRY  TO  CLARKE. 

INDEX=1» 

TANB1  =  B0VRA»(SIN»  (LTD /cos.  (LTD  ) 

TANB2  =  B0VRA»(SIN.  (LT2)/C0S. (LT2)  ) 

C0SB1  =  1./SQRT.  (1.  +  (TANB1*TANBD) 

COSB2=l . /SORT. ( 1 .+ ( TANB2*TANB2 ) ) 

SINB1=TANB1»C0SB1 

SINB2=TANB2»COSB2 

C1=SINB1»SINB2 

Dl=COSBl»COSB2 

DIFL0N=LG2-LG1 

WHENEVER  DI FLON. L  .  0 .  .  I NDEX  =  -1 . 

DIFLON=.ABS.DIFLON 

CDIF=COS.(DIFLON) 

CDIS=C1+D1#CDIF 
SDIS=SORT. ( l.-CDIS*CDIS) 

CA=D1»SIN.(DIFL0N)/SDIS 

CB=CA»CA 

CC=CDIS*( l.-CB) /VK3 
CD=VK4»C1 
CE=VK5»C1 
CF=VK6»CC 

CG1=2»»ATAN» (SORT . ( ( 1 .-CDI S )/(l.+CDIS) ) ) 

CG=CG1/SDI S 

CX=CA*( (CG1#(VK1+CB)+SDIS#( CC+CD)+CG*(CE+CF) ) /VK2) 

DELTAL=CX+DIFLON 

SDELTL=SIN. (DELTAL) 

COELTL=COS. (DELTAL) 

DEN  =  TANB,2»C0SB1-SINB1*CDELTL 


DIR=ATN1.( SDELTL.DEN) 

WHENEVER  01 R.6.P I »DIR=DI R-TPI 
DIR0=DIR7RAD 
DIRD=D1RD*INDEX 
CPH0-C1+D1*CDELTL 
SPH0=SQRT. ( l.-CPH0*CPH0) 

CB0aDl*SDELTL/SPH0 

APH0*2.*ATAN.(SQRT.( ( 1 . -CPHO) / ( 1 .+CPH0 ) } ) 
SB02=1.-CB0*CB0 
C2DEL={2.*C1/SB02 )-CPHO 
C4DEL=(2.*C2DEL*C2DEL)-1. 

SB04=SB02*SB02 
S2PH0=SIN. ( 2.*APH0) 

A0=1 .+E1»SB02-E2*SB04 

B0=E1*SB02-E3*SB04 

C0=E4*SB04 

Drs  =  BRAD*(  AO*APHO+BO*SPHO*C2DEL-CO*S2PHO*C4DEL) 
FUNCTION  RETURN 
END  OF  FUNCTION 


APPENDIX  V 


R  conversion  of  public  land  survey  information 

R  INTO  LATITUDE  AND  LONGITUDE 

R  SUBROUTINES  NEEDED  ARE  DEGRAD*  RADEG.  SPHERE.  AVERAD 
R  W.  R.  TOBLER  /UNIVERSITY  OF  MICHIGAN  /  GEOGRAPHY 
JCOMPILE  MAD. PRINT  OBJECT.PUNCH  OB J ECT . EX ECUT E  TRC 

INTEGER  MER.C1.C2.C3.C4.C5 
INTEGER  COMPAR.N .TWP.RNG .01 .02 .S .PR  INC .N1 
INTEGER  R.S.T 

D'N  SECE(37> .SECN(37) .PMERIDI 36).BLINE(36) .RMER(36) 

VS  DLT(1)=43. 0.43. 0.35. 0.31. 0.36. 0.61. 0.64. 0.34. 0.40. 0.40.0. 

1  42.0.33.0.40.0.34.0.34.0.31.0.42.0.37.0.35.0.34.0.45.0*40.0. 

2  34.0.38.0.60.0.40.0.30.0.30.0.30.0.38.0.40.0.39.0.30.0*45.0. 

3  43.0 

VS  MLT  (11=59.0. 22. 0.1. 0.52. 0.30. 0.49.0. 5  1.0. 38. 0.59. 0.0.0. 

1  30. 0.22. 0.25. 0.59. 0.29. 0.0. 0.25. 0.52. 0.44. 0.15. 0.47. 0.46.0. 

2  7.0.28.0.7.0.0.0.59.0.59.0,26.0*28.0.25.0.6.0.59.0.31.0.0.0 
VS  SLT(l)  =4  4. 0.21. 0.53. 0.32. 0.5. 0.21. 0.50. 0.45. 0.22. 0.50.0. 

1  27.0.38.0.2.0.27.0.32.0.31.0.28.0.54.0.56.0.35.0.13.0.11.0. 

2  20.0.14.0.36.0.7.0.56.0.51.0.3.0.27.0.59.0.23.0.56.0.11.0* 

3  41.0 

VS  DLG(  11=1 04. 0,116. 0.89. 0,90. 0.103. 0.145. 0.1 47. 0.91. 0.84.0. 

1  90.0.90.0,112.0.124.0.86.0,97.0.92.0,84.0.121.0.108.0.106.0. 

2  111. 0,1 11. 0,116. 0.86. 0.1 49. 0.97. 0.91. 0,88. 0.84. 0.89. 0.109.0. 

3  108.0,91.0,122.0,108.0 

VS  MLG( 1)=3. 0.23. 0.14. 0,14. 0,0. 0,18. 0.38. 0,3. 0,48. 0.27.0, 

1  25. 0.18. 0,7. 0.34. 0.14. 0.24. 0,21. 0.54. 0,31. 0,53. 0.39. 0.53.0, 

2  55,0.27.0,21.0,22.0,9.0,1.0,16,0,8.0,56.0,31.0,9.0.44.0,48.0 
VS  SLG(  1)  =16. 0,35. 0,47. 0.41. 0,7. 0.13. 0,26. 0.7. 0,1 1,0, 11.0, 

1  37, 0,19. 0.10. 0,16. 0.49. 0,55. 0,53. 0,47.0. 59, 0.12. 0,33. 0,27.0, 

2  17. 0.21*0.24.0. 8.0. 36. 0.20*0. 38.0, 54.0. 6.0. 59.0.36*0. 34*0. 

3  49.0 
SECE(0  )=0.0 

VS  SECE(l)=2.5,1.5,0.5.-0.5.-1.5,-2.5,-2.5,-1.5,-C.5,0.5, 
11.5.2*5.2.5.1.5.0,5.— 0,5. -1,5. -2*5. -2*5. -1*5. -0.5. 0*5, 1*5. 

22. 5. 2*5. 1,5. 0.5. -0.5. -1.5. -2. 5. -2, 5. -1.5. -0.5. 0.5. 1.5. 2*5 
SECN ( 0  )  =0,0 

VS  SECN(1)=2*5.2.5.2.5.2,5.2.5.2,5.1,5.1.5.1.5.1.5.1*5.1.5. 
10. 5. 0*5. 0*5. 0*5. 0*5. 0.5. -0*5. -0*5. -0.5. -0*5. -0*5, -0*5. -1.5. 

2 -1.5, -1,5, -1.5, -1.5, -1.5. -2. 5, -2. 5, -2. 5, -2. 5, -2. 5, -2. 5 
ra  =  63782064E-04 
MI LE=0. 62136994 
ESQR=6768658E-09 
RAD=174532925E-10 
T'H  INITAL,  FOR  1  =  1. 1, I. G. 35 
DLG (  I  ) =-DLG ( I ) 

EXECUTE  DEGRAD. ( DLT (  I ) , ML T (  I ) , SL T (  I  )  , BL I NE ( I )  ) 

EXECUTE  DEGRAD. ( DLG( I ) ,MLG( I ) ,SLG( I ) .PMERID! I ) ) 

SMLT  =  SIN.( BLINE(  I )  ) 

DUM=SQRT.( 1 .-ESQR*SMLT*SMLT ) 

DUMCUB=DUM,P.3 
DUMMY= ( l.-ESQR )*R1 
RMER ( I ) =DUMMY*MI LE/DUMCUB 
INITAL  CONTINUE 

R'T  CONS.COMPAR 
VS  C0NS=$S3,I1*S 
N  =  0 


START 


N1  =  0 

P'T  SKIP 
VS  SKIP=$1H1*$ 

W'R  C0MPAR.6E.1 

R ' T  LATLON  f Q1 »Q2  »S»T .TWP*R*RN6*C1*C2*C3*C4*C5* 

IDL AT. ML AT. SLAT .DLON .MLON .SLON 
V*  SLATLON=$2Cl .SlO.I2.S3.I2.Cl.S3.I2.Cl.S2.5Cl*S16.F3«0t 
12F2.0.S1.F4«0.2F2.0*$ 

EXECUTE  DEG RAD. ( DLAT .ML AT. SLAT .RLAT ) 

EXECUTE  DEGRAD. ( DLON .MLON . SLON .RLON ) 

O'  E 

R ' T  INDAT.Ql ♦Q2.S.T.TWP.R.RNG.C1.C2.C3.C4.C5 

VSINDAT  =$2C1  .SIO.  12  .S3  .  I  2  .Cl  .S3  .  I  2  .Cl  .S2 . 5C1  *S 

E'  L 

N  =  N+1 

N1=N1+1 

W'R  TWP.E.SNS 

A=T*6. 0-3.0 

O'R  TWP.E.SSS 

A=-(T*6. 0-3.0) 

O'  E 

T'O  ERR 
E'L 

W'R  RNG.E.SES 
B=R*6. 0-3.0 
O'R  RNG.E.SWS 
B=-R*6. 0+3.0 
O'  E 

T'O  ERR 
E'L 

W'R  Ol.E.SNS 

QN=0.25 

O'R  01.E.$S$ 

ON=-0.25 
O'R  Ql.E.S  S 
QN=0.0 
O'  E 

T'O  ERR 
E'L 

W'R  Q2.E.$E$ 

QE=0.25 
O'R  Q2.E.$W$ 

QE=-0.25 
O'R  02.E.S  $ 

QE“0.0 

O'E 

T'O  ERR 
E'L 

W'R  Cl.E.SBS 
W'R  C2.E.SOS 
MER-2 
O'E 
MER=1 
E'L 

O'R  Cl.E.SCS 
W'R  C2.E.$IS 
MER=5 

O'R  C2.E.$OS 
MER  =  6 


O'E 

W'R  C3.E.$IS 
MER  =  3 
O'E 
MER-4 
E'L 
E'  L 

O'R  Cl.E.SFS 
MER  =  7 

O'R  Cl.E.SSS 
MER  =  8 

O'R  Cl. E. SIS 
MER  =  9 

O'R  i..l.E.S^S 

W'R  (C4.E.SAS) .OR. (CS.E.SAS) 

MER-=10 

O'E 

MER=11 

E'L 

O'R  Cl.E.SGS 
MER=12 

O'R  Cl.E.SHS 
W'R  C3.E.SMS 
MER*13 
O'E 

MER«14 

E'L 

O'R  Cl. E. SIS 
MER=15 

O'R  Cl.E.SLS 
MER*16 

O'R  Cl. E. SMS 
W'R  C2.E.SI$ 

MER=>17 

O'E 

MER=18 

E'L 

O'R  Cl.E.SNS 
W'R  C2.E.SA$ 

MER=19 

O'E 

MER=20 

E'L 

O'R  Cl.E.SPS 
MER=21 

O'R  Cl.E.SSS 
W'R  C3.E.SLS 
MER=22 

O'R  Ca.E.SNS 
MER=23 

O'R  ca.E.sws 
MER=25 

0'R(C3.E.SHS).OR.(C4.E.SH5).OR.(C6.E.SH$) 

MER=27 

O'E 

MER=28 
E  '  L 

O'R  C1.E.S2S 
MER=24 


O'R  C1.E.S6$ 

MER=26 

O'R  Cl. E. STS 
MER=29 

O'R  Cl.E.SUS 
W'R  C2.E.S1S 
MER=31 
O'  E 

MER=32 

E'L 

O'R  Cl.E.SWS 
W'R  C3.E.SNS 
MER=35 
O'  E 

MER=34 

W'R  C2.E.SAS.MER=33 
E'L 

O'R  C1.E.S3S 
MER  =  30 
O'  E 

T'O  ERR 
E'L 

A=A+SECN (S)+QN 
A=A/RMER (MER ) 

LATIT=BLINE(MER)+A 
CLAT=COS. ( LATIT  ) 

SMLT=SIN.(LAT1T) 

DUM=SQRT .( l.-ESaR*SMLT*SMLT) 

RPAR=R1*MILE*CLAT/DUM 

B=(B+SECE(S)+QE)/RPAR 

LONGIT=PMERID(MER)+B 

EXECUTE  RAD EG. ( LAT I T .LTD . L TM  .LT  S  ) 

EXECUTE  RADEG. ( LONGIT  »LGD . LGM » LGS ) 

P'T  ONE.Nl 
VS  0NE=S1H 

P'TRI»Q1»Q2.S.T.TWP.R»RNG»C1.C2»C3.C4.C5. 

ILTD.LTM.LTS.LGD.LGM.LGS 

V'SRI=$1H  .2C1.10H  1/4.  SEC  .I2.3H.  r.I2.C1.3H.  R.I2.C1. 

12H.  .S2.5C1.10H.  meridian  / / S 1 . 2 ( F5 . 0 . F 3 . 0 . F3 . 0 . S5 ) . 

219HCALCULATED  LAT/LONG  *S 
W'R  COMPAR.GE.l 
DIFLON=LONGI T-RLON 
DIFLAT=LAT IT-RLAT 

EXECUTE  RADEG. (DIFLAT .DLTD.DLTM.DLTS) 

EXECUTE  RADEG. ( D I FLON .DLGD . DLGM . DLGS ) 

EXECUTE  SPHERE. ( RLAT .RLON .LATIT .LONG  I T  .RHO .ALPHA ) 

EXECUTE  AVERAD. ( RLAT .LAT I T .R1 .ESQR .MRAD  ) 

RHO=RHO*MRAD 

ALPHA=ALPHA/RAD 

P'TR2.DLAT.MLAT.SLAT .DLON .MLON . SLON . DL T D . DL  TM . DLTS . DLGD . 

IDLGM  .DLGS. RHO. ALPHA 

VS  R2  =  51H  ,2  (  F5. 0  .F3 .0  .F3 .0  .  S5  )  .  1 VHOBSERVED  LAT/LONG  /51.2(F 
15.0.F3.0.F3.0.S5) .lOHD I FFERENC E  /S1.F11.4.S6.F9.4.S6. 
222HMILES  AND  DIRECTION  *S 

PUNCH  FORMAT  OUT . N1 . RHO . A LPHA 
V  SOUT  =  $ I 5.S2. Fll .4.S2.F9.4*$ 

E'L 

TRANSFER  TO  START 
PRINT  FORMAT  ONE.N 


ERR 


PRINT  COMMENTS  THIS  OBSERVATION  IS  INCORRECTLY  RECORDEDS 
N1=N1-1 

TRANSFER  TO  START 
E»M 


SCOMPILEMAD.PUNCHOBJECT 

EXTERNAL  FUNCTION  (LLT*ULT*R1 

R  mean  radius  on  ellipsoid 

R  LATITUDES  IN  RADIANS 
ENTRY  TO  AVERAD. 

SMLT=SIN.( ( LLT+ULT)/2. 1 

nuM  =  .';QRT.  (  1  ,-ESQR*SMLT*SMLT  ) 


DUMCUB=DUM*DUM*DUM 

DUMMY= ( 1 .-ESQR 1*R1 
RMER=DUMMY/DUMCUB 


,ESQR»R3) 


RPAR=R1/DUM 
R3=SQRT. (RMER*RPAR) 
FUNCTION  RETURN 
END  OF  FUNCTION 


SCOMP  I  LEMAD  *PUNCHOB  JECT 

EXTERNAL  FUNCTION  ( DEG ♦ M 1 N *  SEC  *  RAD ) 

R  SUBROUTINE  TO  CONVERT  DEGREES  TO  RADIANS 
ENTRY  TO  DEGRAD. 

vector  values  RADI AN= 174532925E-10 
SIGN=RADIAN 

WHENEVER  DEG.L.O.*  SIGN=-RADIAN 

RAD  =  SIGN*(  .ABS.  (DEG)  +  (MIN/60.  )  +  tSEC/36i)0.  )  ) 

FUNCTION  RETURN 

END  OF  FUNCTION 


AVERAD 


DEGRAD 


SCOMPILEMAD.PUNCHOBJECT 

EXTERNAL  FUNCTION  ( R AD .DEG » M I N » SEC ) 

R  CONVERTS  RADIANS  TO  DEGR EES . M I NUTES . 


TMTPnPR  T 


entry  to  radeg. 

VECTOR  VALUES  CONS=206264806E-3 
SEC=.ABS. ( RAD)*CONS 


I=SEC/3600. 
REMAIN=SEC- ( 1*3600. ) 


DEG= 1*1 . 

I =REMAIN/60. 

MI N= 1*1  . 

SEC  =  REMA IN- ( 1*60.  ) 

WHENEVER  R AD .L . 0 . .DEG= -DEG 
FUNCTION  RETURN 
END  OF  FUNCTION 


RADEG 

AND  DECIMAL  SECONDS 


SCOMPILE 


lAD. PRINT  OBJECT. PUNCH  OBJECT 
R  COMPUTES  OBLIQUE  SPHERICAL  ^OO^DINATES 
EXTERNAL  FUNCTION(NLT  *NLG*LAT*  LON  »  RH02  *  GA ) 

VECT0RVALUESP1=314159265E-8 

VECTORVALUESTPI=628318531E-8 

VECTORVALUESPIOVR2=137079633E-8 

VECTORVALUESEPS=0. 0000001 
VECTORVALUESRAD=174532925E-10 


SPHERE 


ENTRY  TO  SPHERE. 

WHENEVER NLT.E. (gO.^RAD) 

GA=LON-NLG 

RH02«PI0VR2-LAT 

OTHERWISE 

WHENEVER (LT.NE.NLT ) .OR. ( LG. NE.NLG ) 

PI=314159265E-8 

TPI=2.*PI 

PIOVR2«P 1/2. 

EPS=0. 0000001 
CNLT=COS. ( NLT ) 

SNLT=SIN. ( NLT ) 

END  OF  CONDITIONAL 
WHENEVER  LON.NE.LONl 
DIF=LON-NLG 
CDIF=COS. ( DIF ) 

SDIF=SIN.(DIF) 

END  OF  CONDITIONAL 
CLT=COS. ( LAT ) 

SLT=SIN. (LAT ) 
0=SLT*SNLT+CLT*CNLT*CDIF 
WHENEVER  O.GE.l. 

RH02=0. 

ORWHENEVER  O.LE.-l. 

RH02=PI 
OTHERWISE 
RH02=ARCC0S. (0) 

END  OF  CONDITIONAL 
NUM=CLT»SDIF 

DEN=CNLT*SLT-SNLT*CLT*CDIF 
WHENEVER . ABS . DEN . L . EPS 
WHENEVER.ABS.NUM.L.EPS 
GA  =  0. 

OTHERWISE 

GA=PI0VR2 

WHENEVER  NUM. L .0. * GA=-GA 
END  OF  CONDITIONAL 
ORWHENEVER. ABS. NUM. L. EPS  •" 

GA  =  0. 

WHENEVER  DEN . L . 0. .GA=P I 

OTHERWISE 

GA  =  ATN1.  (NUM. DEN) 

WHENEVER  GA . G . P I »GA=GA-TP I 
END  OF  CONDITIONAL 
L0N1=L0N 
LT  =  NLT 
LG  =  NLG 

END  OF  CONDITIONAL 
FUNCTION  RETURN 
END  OF  FUNCTION 
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